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Superconductor/ferromagnet/superconductor heterostructures exhibit a so-called long-range 
proximity effect provided some layers of conical magnet Holmium are included in the respective 
interface regions. The Ho layers lead to a spin-flip process at the interface generating equal-spin 
spin-triplet pairing correlations in the ferromagnet. These equal-spin spin-triplet pairing correla¬ 
tions penetrate much further into the heterostructure compared to the spin-singlet and unequal-spin 
spin-triplet correlations which occur in the absence of Ho. Here we present calculations of this ef¬ 
fect based on the spin-dependent microscopic Bogoliubov—de Gennes equations solved within a 
tight-binding model in the clean limit. The influence of the ferromagnet and conical magnet layer 
thickness on the induced equal-spin spin-triplet pairing correlations is obtained and compared to 
available experimental data. It is shown that, in agreement with experiment, a critical minimum 
thickness of conical magnet layers has to be present in order to observe a sizeable amount of equal- 
spin spin-triplet pairing correlations. 

PACS numbers: Valid PACS appear here 


I. INTRODUCTION 

The proximity effect arises because superconducting 
pairing correlations are able to penetrate the interface be¬ 
tween a superconductor (SC) and a nonmagnetic metal. 
The correlations decay in the adjacent metal with a char¬ 
acteristic length scale. Replacing the nonmagnetic metal 
by a ferromagnet (FM) the pairing correlations become 
oscillatory and this characteristic length scale is also con¬ 
siderably reduced. The penetration depth of the spin- 
singlet pairing correlation into the ferromagnetic layer 
depends on the exchange interaction, which also leads to 
FFLO-like oscillationsi*^ visible in the ferromagnetic re¬ 
gion. The SC/FM interface also generates unequal-spin 
spin-triplet pairing correlations which have the same os¬ 
cillating and decaying pattern as for the spin-singiet cor¬ 
relations. However, according to a theoretical predic¬ 
tion by Bergeret et al^ it should be possible to also cre¬ 
ate equal-spin spin-triplet pairing correlations, provided 
the interface allows for some kind of spin-flip process. 
In principle, these equal-spin spin-triplet pairing corre¬ 
lations should be compatible with a ferromagnetic ex¬ 
change interaction and allow for much larger penetration 
depth in the ferromagnetic region of the heterostructure 
compared to the spin-singlet proximity effect. This so- 
called long-range proximity effect has stimulated a lot of 
experimental and theoretical work, and has been sum¬ 
marised in reviews by Buzdini and Bergeret et al^ 

From the experimental point of view several possible 
sources of spin-flip processes have been identified and 
have been realised in heterostructure setups. It should 
be noted that experimental evidence for the long-range 
proximity effect mostly stems from observations of su¬ 
per currents in SC/FM/SC Josephson junctions with FM 
length scale incompatible with the spin-singlet proximity 
effect. Experimentally the heterostructures which have 


been already realised include half-metallic metals)^^— in¬ 
troducing a magnetic inhomogeneity at the interface,— 
noncollinear magnetic interfaces'^— and helicali^ or 
conical magnetic structures^^^— to name but a few. 

On the theoretical side, several basically very different 
approaches have been applied to particular heterostruc¬ 
tures to investigate equal-spin spin-triplet pairing cor¬ 
relations and how they are affected by specific spin-flip 
mediating interfaces. In line with the experimental ob¬ 
servations various interfaces which have been studied in¬ 
clude half-metallic metalsj ^^d^ inhomogeneous magneti- 
sations;^*^^ noncollinear magnetisations^^ and helical 
(or conical) magnets— Green’s function techniques 
based on solutions of the Eilenberger— and the Usadel 
equation have been reported )^d^i^°'^^i^^ as well as self- 
consistent solutions of the Bogoliubov—de Gennes (BdG) 
equations for suitable tight-binding models— 

Here we present results based on a model heterostruc¬ 
ture consisting of a s-wave SG/FM/SG junction with ad¬ 
ditional conical magnet (CM) layers introduced at the 
interfaces. The relevant tight-binding model is solved in 
the microscopic spin-dependent BdG equations and solu¬ 
tions are iterated to self-consistency. From the respective 
eigenfunctions we obtain the different spin-triplet pair¬ 
ing correlations and discuss the influence of CM and FM 
layer thickness. We find that, in agreement with experi¬ 
mental observations, a critical minimal thickness of CM 
layers has to be present to observe the long-range prox¬ 
imity effect. 

The paper is organised as follows. Sec. HIl provides the¬ 
oretical background, including the BdG equations and 
heterostructure setup in Sec. Ill A[ and the spin-triplet 
pairing correlations in Sec. IHBl respectively. Sec. Hill 
is devoted to results, where the influence of FM (CM) 
thickness on spin-triplet pairing correlations is discussed 
in detail in Sec. IHI Al (Sec. IHIBI) . A summary and con¬ 
cluding remarks are given in Sec. IIVI 
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II. THEORETICAL BACKGROUND 


A. Method, heterostructure setup and computational details 


The results presented here are based on self-consistent solutions of the microscopic BdG equations in the clean 
limit. For the spin-dependent case and incorporating the vector components of a general exchange field h the BdG 
equations rea d^^i^^i^^ 
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with e„, and Una- and Vna denoting the eigenvalues, and 
quasiparticle and quasihole amplitudes for spin ct, respec¬ 
tively. The tight-binding Hamiltonian 'Ho can be simpli¬ 
fied according to ^^i^° 

Ho = —t 'y ] ^cl^Cn+l + Cn^ + ^ ^ ~ M) (2) 

n n 

At multilayer index n we have the electronic creation 
(cjj) and destruction operators (c„), whereas the next- 
nearest neighbour hopping parameter and the chemical 
potential (Fermi energy) are chosen to be t = 1 and ^ = 
0, respectively. 

According to Balian and Werthamer— >2^ the general 
form of the pairing matrix in Eq. © can be rewritten as 


where we are summing only over positive eigenvalues e„ 
and /(£n) denotes the Fermi distribution function evalu¬ 
ated as a step function for zero temperature. Setting up 
the multilayer structure as shown in Fig. [TJa) the effec¬ 
tive superconducting coupling parameter g{r) equals 1 in 
the nsc = 250 layers of spin-singlet s-wave superconduc¬ 
tor to the left and right of the heterostructural setup and 
vanishes elsewhere. 

The thickness of the CM layer ncM is varied from 0 
to 25 layers, and we include up to npM = 500 layers 
of FM in the middle of the heterostructure. The vector 
components of the conical exchange field are determined 

by 2 i >22 
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with • • • indicating a 2 x 2 matrix, respectively. Making 
use of the Pauli matrices d, the superconducting order 
parameter is described by a singlet (scalar) part A and 
a triplet (vector) part d, respectively. Defining A to 
contain only the triplet d-vector components yields 


A = (dd) ia2 = 


-dx + idy dz 


idy 


( 4 ) 


h = /iQ s cos ay + sin a 
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( 7 ) 

with do = 0.1 fixing the exchange field in the CM layers, 
and a = 1 being the lattice constant. The opening angle 
a (measured from +y towards +z) and the turning angle 
(3 (measured from +z towards -l-x) are fixed to the values 
of Ho, namely a = 80 ° and /3 = 30 °. 


Considering this general form of this triplet pairing ma¬ 
trix the product AA^ can be written as 

AA’^ = |dpdo -I- J (d X d*) d, (5) 

with |d| describing the gap function and d x d* being 
a measure for the Cooper pair spin magnetic moment, 
respectively. 

However, here we are interested in the case of s-wave 
SC only. Therefore the pairing potential is restricted to a 
scalar quantity A fulfilling the self-consistency condition 

- f[£n)] 

( 6 ) 

+ Mnt(rX^(r)/(e„)), 


*^SC ^CM ^FM 



FIG. 1. (a) Heterostructure setup (from left to right): spin- 
singlet s-wave superconductor (nsc layers), conical magnet 
(ncM layers), ferromagnetic metal (npM layers), conical mag¬ 
net (ncM layers), and spin-singlet s-wave superconductor 
(nsc layers), (b) Opening and turning angles a and /3 of 
the conical magnet. According to Eq. 0 Q is measured from 
-\-y towards -\-z, whereas /3 is measured from +z towards -\-x. 
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B. (Triplet) Pairing correlations 


The superconducting pairing correlation between spins 
a and /3 can generally be evaluated as on-site average for 
times t = T and t' = 0 as 


/a/3(r, r, 0) = 0)), (8) 

I 


+ fir(y^T)) = ^ 
h{y,r) = ^{fniy^'r) - fuiv^T)) = ^ 


depending on position y and time parameter r (fixed to 
r = 10 throughout this work). The r dependence is 
governed via C„ (r) given by 

Cn(T) = cos(£r„T) - zsin(£r„T)(l - 2/(e„)). (10) 

The different pairing correlations can be rewritten simi¬ 
larly to the pairing matrix in Eq. 


f fn /ti 
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— fo + fz A + i/y / 


( 11 ) 

In analogy to Eq. o a restriction to the vector com¬ 
ponents of the f-vector yields the triplet pair function 
matrix^^ 


Note, that we introduced an additional factor i compared 
to-^ to be consistent with the definition of A in Eq. Q. 
Instead of Eq. m we now have 

/A = |fpdo + *(fxr)d, (13) 

with the two analogues of |d| and d x d*, being |f| and 
f X f*, being conveniently expressed in terms of the f- 
vector components 


A — 2 (“At + At) 

A = -^(At + At) (14) 

A = 2 (At + At) ■ 


with d/o-(r, r) being the many-body field operator for spin 
a at time r. The time-dependence is introduced through 
the Heisenberg equation of motion. A pairing correlation 
evaluated using Eq. m is local in space, leading to van¬ 
ishing triplet contributions for t = 0 in accordance with 
the Pauli principlei^ The nonvanishing contributions for 
finite times r are an example of odd-frequency triplet 
pairingi^ Substituting the field operators valid for our 
setup and phase convention the spin-dependent triplet 
pairing correlations read 


iurit(.y)Vni(.y) + Uni{y)Vn^{y))Cn{T) 
n 
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III. RESULTS AND DISCUSSION 


A. Influence of ferromagnet thickness ufm 


The self-consistent calculation as described in Sec. El 
yields the spin-triplet pairing correlations A and A as 
defined in Eq. (|9l). For an exemplary heterostructure 
setup with npM = 100 layers and one full turn of the 
CM to either side, the upper (lower) left panels of Fig. [2] 
show results on the real (green) and imaginary (orange) 
parts of the spin-triplet pairing correlations A (fi), re¬ 
spectively. The unequal-spin spin-triplet pairing correla¬ 
tion A clearly shows the expected oscillatory behaviour 
in the EM region of the heterostructure. These oscilla¬ 
tions change in the CM region and decay in the adjacent 
SC layers. In contrast, the equal-spin spin-triplet pairing 
correlations A show pronounced features only inside the 
CM region and decay into both the EM and the adjacent 
SC layers. However, keeping in mind that A IIS'S contri¬ 
butions from At snd Ai (see Eq. ([H])), the upper (lower) 
middle panels of Fig. [5] show the real (green) and imagi¬ 
nary (orange) parts of At (At) separately, clearly show¬ 
ing a nonvanishing contribution inside the EM region for 
the single spin-channels. A more detailed analysis of the 
influence of CM orientation in the heterostructure on the 
different spin-triplet pairing correlations can be found in 
an earlier worki^^ Recalling the pair function matrix as 
of Eq. (fT51l the upper and lower right panels of Fig. [2] 
show the contributions to the magnitudes of the f-vector 
and f X f*, respectively. The magnitude of the f-vector 
resembles the oscillatory patterns already known from 
A (upper left panel). The magnitude of f x f*, associ¬ 
ated with the spin magnetic moment of the Cooper pairs, 
shows leakage from the CM layer into the SC layer. How¬ 
ever, this leakage shows a strong decay and vanishes after 
a few layers. On the contrary, there is no such leaking 
from the CM into the EM region of the heterostructure. 
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FIG. 2. (colour online) Real (green) and imaginary (orange) parts of the spin-triplet pairing correlations /o (upper left panel) 
and fi (lower left panel) according to Eq. ([9]). For better analysis the upper (lower) middle panels depict the real (green) and 
imaginary (orange) parts of /-|-|- (fn) as components of fi shown in the lower left panel. Upper and lower right panels show 
the magnitude of the f-vector and the magnitude of f x f* as introduced in Eq. m, respectively. All data is shown depending 
on the multilayer index n, with n = 0 lying in the centre of the ferromagnetic layer. The vertical dashed lines indicate the 
FM/CM and CM/SC interfaces, respectively. 




Keeping the number of CM layers fixed to one full turn 
we now vary the number of FM layers from npM = 0 to 
npM = 500 layers. The magnitudes of /o and fi depend¬ 
ing on npM are depicted in the upper and lower panels 
of Fig. O Similarly to the results presented in Fig. [21 
the upper panels in Fig. [31 depicting the magnitudes of 
/o, show a pronounced oscillatory behaviour within the 
FM region up to npM = 500 layers. There are also much 
weaker oscillations visible inside the SC region depend¬ 
ing on the number of FM layers considered. As expected 
from the lower left panels of Fig. [2] there is no contri¬ 
bution to the magnitude of f\ in the SC or FM region 
of the heterostructure. The only sharp signals originate 
from the CM regions and are of similar strength as /q. 
The strengths of both, /o and /i, are influenced by the 
number of FM layers npM, especially visible for small 
values of npM, where oscillations are seen in /o and fi. 
Recalling the pair function matrix as of Eq. (ITSl) the mag¬ 
nitudes of the f-vector and f x f * are shown in the upper 
and lower panels of Fig. 01 respectively. As can be seen 
from the upper panels the magnitude of the f-vector is 
quite similar to the magnitude of /o already presented in 
Fig. 121 Keeping in mind the appearance of the magni¬ 
tude of f X f* shown in the lower right panel of Fig. |21 
there is no influence of the number of FM layers npM on 
this equal-spin pairing correlation. For all npM consid¬ 
ered here, the strongest contributions to the magnitude 
of f X f * stem from the CM layers only, with smaller and 
fast decaying contributions in the adjacent SC layers. No 
contributions are visible inside the FM layers, as already 
noted in the discussion of Fig.|21 To summarise this sec¬ 
tion, the number of FM layers npM in the heterostruc¬ 
ture has no influence on the magnitudes of fi and f x f *, 
whereas the number of maxima in the magnitudes of /o 


and the f-vector are increased with increasing number of 
FM layers npM- 


B. Influence of conical magnet thickness ncM 

Understanding the behaviour of FM layers on the 
various spin-triplet pairing correlations as discussed in 
Sec, umi we are now interested in the influence of the 
number of CM layers ncM- For the following calculations 
we fix the number of FM layers to npM = 100 layers and 
vary ncM from 0 to 25 layers, respectively. The mag¬ 
nitudes of /o and /i depending on ncM are depicted in 
the upper and lower panels of Fig. jS] Keeping in mind 
that npM is fixed to 100 the upper panels of Fig. [5] show 
the already familiar oscillatory patterns in the magnitude 
of /q. In contrary to previous results, the magnitude of 
/o is affected by the number of CM layers ncM, as can 
be clearly seen in the intensity variations in the FM re¬ 
gion shown in the upper right panel of Fig. [S] As the 
number of CM layers ncM increase, the magnitudes of 
both, /o and /i, show a strongly oscillatory behaviour 
in the CM region of the heterostructure. Minima ap¬ 
pear at multiples of half-integer turns of the full conical 
magnetic structure, whereas the maxima lie inbetween. 
Also, there is a phase shift present between the maxima 
appearing in the magnitudes of /o (upper left panel of 
Fig. |S]) and fi (lower left panel of Fig. [5]). In addition, 
the magnitude of fi show nonvanishing contributions in 
the FM region of the heterostructure, showing the same 
oscillatory dependence as in the CM regions, respectively. 
Apparently, at least a quarter of a full conical magnetic 
structure has to be present in the heterostructure to ob¬ 
tain a sizable effect on the magnitude of fi (lower panels 
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FIG. 3. (colour online) Magnitudes of /o and /i depending on number of FM layers riFM varied from tifm = 0 to ufm = 500 
layers, respectively. Full data sets are shown in the left panels, whereas the right panels depict a top view of the data. 




FIG. 4. (colour online) Magnitudes of the f-vector and f x f* as introduced in Eq. (1131) depending on the number of FM layers 
riFM varied from u-fm = 0 to ufm = 500 layers, respectively. Full data sets are shown in the left panels, whereas the right 
panels depict a top view of the data. 


of Fig. [5]). The respective magnitudes of the f-vector and 
f X f* are shown in the upper and lower panels of Fig. [51 
respectively. With the magnitude of fi now showing con¬ 
siderable contributions, the difference between the mag¬ 
nitudes of /o (upper panels of Fig. [5]) and the f-vector 


(upper panels of Fig. |6l) is larger compared to Sec. IIII Al 
The magnitudes of /i are added to those of /o to yield 
the magnitudes of the f-vector, and are thus responsible 
for the small intensity variations in the FM region (upper 
right panel of Fig. El). With the slight phase shift present 
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FIG. 5. (colour online) Magnitudes of fo and /i depending on number of CM layers ncM varied from ncM = 0 to ncM = 25 
layers, respectively. Full data sets are shown in the middle panels, whereas in the left and right panels the same data is shown 
when viewed from the yz plane (right CM region only) and as a top view, respectively. 


between the maxima appearing in the magnitudes of fo 
and /i the features visible in the magnitude of the f- 
vector (lower left panel of Fig. |6]) are less sharp. From 
the lower panels of Fig. [HI depicting the magnitudes of 
f X f*, we observe a strong influence of the number of 
CM layers ncM on the intensity. This again shows os¬ 
cillatory behaviour related to the number of turns of the 
conical magnetic structure. However, in this case there is 
no increase in the penetration of f x f * into the FM or SC 
region. Summarising this section, the number of CM lay¬ 
ers ncM has a more pronounced effect on the spin-triplet 
pairing correlations compared to the number of FM lay¬ 
ers npM as discussed in Sec. IIII Al The oscillations in the 
magnitudes of /i shown in the FM and CM regions of 
the heterostructure are also visible in the magnitude of 
the f-vector. In addition, the first sizable contributions 
to the magnitudes of /i and f x f* require at least a 
number of CM layers corresponding to a quarter of a full 
conical magnetic turn. 

IV. SUMMARY AND OUTLOOK 

In summary, we presented results on the spin-triplet 
pairing correlations in a SC/FM/SC heterostructure in¬ 
cluding Ho as a spin-flip source at the interfaces. The 


calculations were based on self-consistent solutions of the 
microscopic BdG equations in the clean limit. In addi¬ 
tion to the spin-triplet pairing correlations fo and /i, the 
magnitudes of the f-vector and f x f * have also been cal¬ 
culated. These allow for a deeper understanding of the 
influence of varying thicknesses of the FM and CM lay¬ 
ers on the total and equal-spin spin-triplet correlations, 
respectively. The leaking of spin magnetic moment of 
the Cooper pairs into the SC region of the heterostruc¬ 
ture is limited and decays rapidly after only a few layers 
of SC. This behaviour is not influenced by the thickness 
of FM or CM layers. On the contrary, sizeable contri¬ 
butions to the magnitude of fi require a CM thickness 
corresponding to at least a quarter of the full magnetic 
cone structure. This is in agreement with experimental 
observationsi^ where similar heterostructures have been 
investigated. Thereby, spin-triplet super currents exhibit 
peak values at CM thicknesses corresponding to roughly 
half a turn of the conical magnetic structure. 
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FIG. 6. (colour online) Magnitudes of the f-vector and f x f* as introduced in Eq. (I13II depending on the number of CM layers 
ncM varied from ncM = 0 to ncM = 25 layers, respectively. Full data sets are shown in the middle panels, whereas in the left 
and right panels the same data is shown when viewed from the yz plane (right CM region only) and as a top view, respectively. 
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